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In  this  paper,  we  investigate  an  adaptive  energy  management  and  power  splitting  system  for  a  fuel  cell 
hybrid  electric  vehicle.  The  battery  pack  is  the  main  power  source  whereas  the  fuel  cell  is  considered  as  a 
range  extender  that  cannot  sustain  alone  the  vehicle  traction  power.  In  addition,  the  fuel  cell  contributes 
to  reduce  the  battery  pack  degradation  by  limiting  its  depth-of-discharge  (DoD).  This  energy  manage¬ 
ment  system  is  based  on  a  two  layer  architecture  in  which  the  upper  layer  computes  the  anticipated  end- 
of-trip  DoD  using  online  mass  estimation.  The  lower  layer  is  designed  to  split  the  driver  power  demand 
by  minimizing  a  cost  function  which  includes  the  hydrogen/electricity  cost  ratio.  Therefore,  the  best 
trade-off  between  reducing  battery  pack  degradation  and  using  cost  effective  energy  is  provided. 
Furthermore,  the  system  allows  the  fuel  cell  to  operate  at  its  maximum  efficiency.  Comparative  study 
results  indicate  that  using  online  mass  estimation  improves  the  overall  fuel  consumption  efficiency 
whilst  contributing  at  the  same  time  to  DoD  reduction. 

©  2013  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Terrestrian  vehicle  powertrain  electrification  is  a  key  step  to¬ 
wards  the  reduction  of  the  greenhouse  gas  emission  sources  [1], 
However,  the  use  of  batteries  as  a  power  source  for  vehicles  raises 
several  challenges  that  must  be  addressed  to  make  it  competitive: 
limited  operating  range,  long  charging  time,  limited  lifespan,  etc. 
Vehicle  source  hybridization  has  been  proposed  as  a  practical  way 
to  increase  fuel  usage  efficiency  and  to  extend  operating  range  [2], 
Thus,  plugin  hybrid  vehicles  (PHEV)  with  gasoline  internal- 
combustion  engines  as  range  extenders  have  been  proposed  by 
car  makers  and  several  of  such  vehicles  are  being  deployed  around 
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the  word:  Chevrolet  Volt  of  General  Motor,  Outlander  PHEV  of 
Mitsubishi,  etc.  However,  it  is  noteworthy  that  gasoline  used  in 
these  vehicles  contributes  to  the  greenhouse  gas  emission  [3,4], 

To  completely  remove  any  fossil  energy  on  the  vehicle,  a  clean- 
energy  vector  such  as  hydrogen  attracts  a  lot  of  attention  in 
research  community  [5j.  Indeed,  hydrogen  has  one  of  the  highest 
energy  densities  per  weight  and  can  be  used  in  fuel  cell  to  produce 
electricity  or  burned  in  an  internal  combustion  engine  (ICE)  to 
directly  generate  propulsion  torque  [6],  Some  of  PHEVs  with 
hydrogen  as  energy  vector  have  been  proposed  in  the  literature  [7— 
9],  Fuel  cell  PHEV  (FC-PHEV)  has  received  more  attention  than 
hydrogen-based  ICE. 

The  vehicle  global  fuel  saving  is  thoroughly  related  to  the  power 
flow  control  between  each  energy  source  and  the  powertrain 
[10,11],  Several  energy  management  and  power  sharing  control 
systems  for  hybrid  electric  vehicles  (HEV)  as  well  as  PHEV  have 
been  reported  throughout  the  literature:  rule-based  methods 
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[8,12],  optimization  methods  [13—15],  stochastic  approaches  and 
adaptive  approaches  [16,17],  These  methods  belong  to  the  blended 
mode  energy  management  systems.  Most  of  these  works  assumed 
that  the  power  demand  from  the  vehicle  driver  is  known.  However, 
the  driver  often  provides  a  control  signal  using  the  vehicle  accel¬ 
eration  pedal.  One  way  to  determine  the  corresponding  power 
demand  is  to  use  the  vehicle  longitudinal  dynamics  which  con¬ 
siders  that  the  mass  and  other  physical  parameters  (rolling  resis¬ 
tance  coefficient,  road  grade,  air  density,  vehicle  frontal  area,  and 
aerodynamic  drag  coefficient)  as  known  beforehand.  Whilst  some 
of  these  parameters  can  be  assumed  to  be  constant  values,  the 
vehicle  mass  as  well  as  the  road  grade  can  change  during  the  trip.  In 
addition,  most  of  the  reported  work  assumed  that  the  on-board 
energy  sources  have  similar  power  range,  and  the  energy  cost  is 
not  directly  taken  into  account.  Having  similar  operating  power 
range  simplifies  the  design  of  the  energy  management  system. 
Indeed,  the  well-known  charge-depleting  charge-sustaining 
method  can  be  easily  used:  after  the  battery  pack  has  been 
depleted,  the  secondary  energy  source  can  sustain  the  power  de¬ 
mand  so  that  at  the  end  of  the  trip,  the  battery  pack  energy  remains 
close  to  a  given  low  value.  In  addition,  in  blended  mode,  the  power 
demand  needs  to  be  estimated  in  order  to  properly  use  the  sec¬ 
ondary  energy  source  to  supplement  the  battery  pack. 

This  paper  investigated  the  energy  management  and  the  real¬ 
time  power  splitting  problems  of  a  serial  topology  FC-PHEV  in 
which  the  battery  pack  is  the  primary  and  the  most  powerful  en¬ 
ergy  source  and  the  fuel  cell  is  considered  as  a  range  extender 
(fewer  powerful  energy  source).  Therefore,  the  charge-depleting 
charge-sustaining  method  cannot  be  used  since  the  fuel  cell  is 
unable  to  sustain  the  power  demand.  So,  to  maintain  battery  pack 
energy  close  to  a  prescribed  minimum  value  at  the  end  of  a  trip,  the 
fuel  cell  operating  sequence  must  be  carefully  designed.  In  addi¬ 
tion,  we  assume  that  when  the  vehicle  is  stopped  during  a  trip,  the 
mass  may  change  (change  in  the  number  of  passengers  in  the 
vehicle,  garbage  trucks,  delivery  trucks  and  buses,  etc.).  Further¬ 
more,  it  can  be  used  for  small  vehicles  that  cannot  support  the 
weight  of  a  great  fuel  cell  (mass-power  ratio). 

The  mass  estimation  in  the  context  of  hybrid  electric  vehicle 
energy  management  has  not  been  fully  addressed  in  the  literature. 
It  becomes  a  key  step  towards  a  good  energy  planning  system,  and 
its  online  estimate  is  challenging.  Several  mass  and  grade  estima¬ 
tion  methods  reported  throughout  the  literature  can  be  classified  in 
two  categories  [18]:  event-based  and  averaging.  Each  of  these 
categories  could  be  further  divided  into  two  sub-category:  simul¬ 
taneous  estimation  of  mass  and  grade,  single  estimation  of  mass 
using  grade  measurements.  For  these  estimations,  different  vehicle 
dynamic  models  were  used  [18]:  suspension  dynamics  [19],  yaw 
dynamics  [20],  drive-train  dynamics  and  longitudinal  dynamics 
[21-23], 

Event-based  methods  seek  for  driving  conditions  that  provide 
sufficient  excitation  for  mass  or  grade  estimation:  model  predictive 
approach  [24],  use  of  grade  estimation  with  GPS  (Global  Positioning 
System)  to  derive  mass  estimation  [25]  and  supervisory  control 
approaches  [18], 

Instead  of  seeking  for  events,  the  averaging  method  continu¬ 
ously  monitors  the  vehicle  dynamics  in  order  to  directly  estimate 
the  mass  and  grade  online,  using  recursive  least  squares  (RLS) 
[21,23,26,27]  and  Kalman’s  filter.  Others  reported  methods  can  be 
seen  in  Refs.  [22,28:. 

The  proposed  energy  management  system  is  based  on  a  two- 
layer  architecture  in  which  the  supervisor  (high  level  of  the  archi¬ 
tecture)  is  responsible  for  providing  the  global  energy  consumption 
profile  when  the  vehicle  mass  is  not  constant  and  the  anticipated 
battery  depth-of-discharge  (DoD).  Given  the  global  consumption 
profile,  the  lower  level  provides  a  locally  optimal  power  sharing 


between  the  fuel  cell  and  the  battery  pack  by  minimizing  a  cost 
function  which  includes  the  hydrogen  and  electricity  price  ratio. 

The  rest  of  the  paper  is  organized  into  seven  sections.  The  FC- 
PHEV  description  and  its  longitudinal  model  are  presented  in 
Sections  2  and  3,  respectively.  Based  on  this  model,  a  new  adaptive 
energy  management  system  is  described  and  discussed  in  Section 
4.  Section  5  is  related  to  the  design  of  the  lower  layer  (power 
splitting  method).  The  adaptive  energy  planning  and  the  power 
splitting  method  validations  are  provided  in  Sections  6  and  7, 
respectively.  Finally,  the  conclusion  is  presented  in  Section  8. 

2.  FC-PHEV  description 

Fig.  1  represents  a  serial  topology  of  a  FC-PHEV  where  the  fuel 
cell  is  not  required  to  follow  precisely  the  dynamics  of  the  driver 
power  demand  if  the  battery  pack  is  well  sized.  In  addition,  with 
this  topology,  the  fuel  cell  can  be  optimally  set  to  provide  power. 

We  consider  that  the  sub-system  comprising  the  uni-directional 
DC— DC  converter  and  the  fuel  cell  represents  the  fuel  cell  power 
source.  The  driver  provides  two  different  signals:  the  brake  and 
acceleration  commands.  Given  these  signals,  the  Power  Demand 
Module  generates  the  corresponding  mechanical  power  request  Pu 
which  is  further  sent  to  Energy  Management  System  (EMS).  Pu  is 
interpreted  as  the  desired  mechanical  power  during  a  trip.  In 
addition,  the  EMS  receives  the  stored  hydrogen  energy  Efc,  the 
available  energy  in  the  battery  pack  £b,  the  different  power  source 
efficiency  maps  (not  represented  in  Fig.  1 )  and  the  hydrogen/elec¬ 
tricity  cost  ratio  (the  cost  of  1  kW  h  of  hydrogen  divided  by  the  cost 
of  1  kW  h  of  electricity). 

Taking  into  account  all  this  information  and  considering  that  the 
vehicle  mass  M  can  vary  during  the  trip,  the  role  of  the  EMS  is  to 
find  the  most  appropriate  power  splitting  of  Pu  by  providing  the 
fuel  cell  optimal  reference  power  command  P£c  and  the  Propulsion 
System  reference  power  command  f*,.  The  fuel  cell  DC— DC  con¬ 
verter  is  responsible  for  providing  Pfc  that  corresponds  to  Pfc.  On  the 
other  hand,  the  Propulsion  System  controller  (not  represented  in 
Fig.  1)  generates  the  vehicle  mechanical  power  Pm  using  the 
reference  command  P^.  Pt,  represents  the  battery  pack  power:  a 
positive  value  indicates  that  the  battery  pack  is  providing  electrical 
power  whilst  a  negative  value  indicates  that  the  battery  pack  is 
being  recharged. 

In  this  paper,  we  consider  the  two-layer  architecture  proposed 
in  Ref.  [29]  in  which  the  upper  layer  is  responsible  of  the  globally 
optimal  energy  profile  during  the  trip  whilst  the  lower  layer  role  is 
to  split  the  power  demand  so  that  the  vehicle  energy  consumption 
follows  this  globally  optimal  profile. 

3.  FC-PHEV  longitudinal  model 

Assume  that  the  wind  speed  relative  to  the  ground  is  negligible. 
The  FC-PHEV  longitudinal  dynamics  is  represented  by  Equation  ( 1 ) 
[14]  and  the  corresponding  mechanical  power  is  given  by  (2). 

Fm(k)=Mv(k)+^p3C<iAv2{k)+Mgsm(d(k))+Mggcos(d{k))  (1) 

where  k,  Fm,  M,  g,  p3,  p,  Q  and  A  represent  respectively  the  sampling 
index,  the  resultant  mechanical  force,  the  vehicle’s  total  mass,  the 
gravity  constant,  the  air  density,  the  rolling  resistance  coefficient, 
the  air  drag  coefficient  and  the  active  frontal  area;  d,  i>  and  v 
represent  respectively  the  road  grade,  the  longitudinal  acceleration 
and  the  vehicle  speed.  0  <k<N  where  N  is  the  number  of  samples 
in  a  given  driving  cycle. 

Pm  (k)  =  Fm(k)v(k) 


(2) 
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Hydrogen 
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Fig.  1.  Fuel  cell  plugin  hybrid  electric  vehicle:  the  arrows  with  the  solid  lines  represent  the  energy  flows;  the  arrows  with  the  dashed  lines  indicate  the  control  and  measured 


The  power  balance  at  the  electric  node  (see  Fig.  1)  is  given  by 
Equation  (3). 


be  determined  through  any  classical  system  identification 
methods. 


Pb(k)  =  Pe(k)  —  Pfc(k) 


(3) 


4.  Energy  management  system 


where  Pe  is  the  electric  power  provided  to  the  powertrain. 

Knowing  the  powertrain  efficiency  map  7jm  and  the  mechanical 
power  at  the  vehicle  wheels,  the  expression  of  Pe  yields: 

Pe(k )  =  ^(k)Pm(k)  (4) 

The  variation  of  battery  pack  energy  Eb  is  represented  by  the 
following  equation: 

Eb(k+  1)  =  Eb(k)  -  Vb\k)Pb(k)AT  (5) 

where  AT  and  i]b  represent,  the  sample  time  period  and  the  battery 
pack  efficiency,  respectively. 

Assume  that  the  fuel  cell  parameters  are  adequately  controlled. 
The  energy  balance  in  the  hydrogen  tank  [30-32]  is  given  by: 

Efc(k  +  1)  =  Efc(k)  -  KH2mH2(k)AT  (6) 

where  EfC,  Kh2  and  jtih2  represent  the  hydrogen  energy  in  the 
storage,  the  hydrogen  high  heating  value  [33,34]  and  the  hydrogen 
flow  rate,  respectively.  The  fuel  cell  power  system  efficiency  r)fc  is 
defined  as: 


Vfc  (k)  = 


Pfc(k) 

Khj'Ph  2(k) 


(7) 


To  take  into  account  the  fuel  cell  dynamics,  the  following 
equation  previously  proposed  in  Refs.  [29,30,35]  is  used: 


Pfc(k  + 1)  =  A(k)PfC(k)  +  B(k)PfC(k) 


(8) 


where  0  <  A(k)<l,  B(k)  and  Pj?c(k)  represent  respectively  the 
Fuel  Cell  dynamic  constant,  the  input  coefficient  and  the 
input  power  command  (see  Fig.  1 ).  The  sequences  of  A  and  B  can 


The  FC-PHEV  has  two  energy  sources:  the  battery  pack  is  acting 
as  the  primary  source  and  the  fuel  cell  is  the  secondary  source  that 
cannot  sustain  alone  the  power  demand.  For  long  range  operation, 
the  battery  pack  which  is  initially  fully  charged,  is  unable  to  provide 
the  total  required  energy.  Instead,  the  fuel  cell  supplements  the 
power  provided  by  the  primary  source.  However,  the  hydrogen  cost 
is  generally  higher  than  the  electricity  from  the  grid  and  this  cost 
ratio  must  be  taken  into  account  when  designing  the  Energy 
Management  System. 

To  allow  the  battery  pack  to  not  being  depleted  below  a  given 
minimum  energy  threshold,  the  EMS  must  determine  the  appro¬ 
priate  time  to  start  the  fuel  cell.  Indeed,  if  the  fuel  cell  is  started 
late  during  the  trip,  the  battery  pack  energy  may  drop  under  this 
minimum  threshold  at  the  end  of  the  trip.  On  the  other  hand,  as 
the  battery  pack  is  assumed  to  be  initially  completely  charged, 
starting  too  earlier  the  fuel  cell  may  reduce  its  capacity  to 
adequately  handle  the  regenerative  power  source  during  the  trip. 
Furthermore,  the  EMS  must  take  into  account  the  fact  that  the 
vehicle  and  load  mass  can  vary  during  the  trip. 

The  upper  layer  which  is  part  of  the  EMS  architecture,  provides 
the  global  and  optimal  energy  consumption  profile.  The  energy 
profile  will  therefore  be  used  to  estimate  the  anticipated  battery 
pack  depth-of-discharge.  In  the  following  section,  this  layer  prob¬ 
lem  formulation  and  design  are  proposed. 

4.1.  Global  energy  planning 

The  first  task  of  the  upper  layer  within  the  EMS  is  to  estimate  the 
global  energy  consumption  for  a  given  trip  by  considering  only  the 
prescribed  speed  limit  and  the  position  of  traffic  lights  and  stop 
signs.  In  addition,  it  must  take  into  account  the  mass  variation. 
Speed  limit  and  position  of  traffic  lights  and  stop  signs  are 
commonly  available  for  each  urban  road  and  highway  segment,  and 
are  independent  of  the  vehicle  driver’s  behavior  [36].  This  two 
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information  and  the  road  grade  can  be  also  retrieved  from  any 
Geographic  Information  System  (GIS)  device  or  any  appropriate 
traffic  flow  database  system. 

The  speed  limit  profiles  for  all  road-trip  segments  may  have 
abrupt  changes  in  the  speed  values  which  may  produce  high  ac¬ 
celerations  if  these  profiles  are  used  as  driving  cycle  for  a  vehicle. 
These  high  accelerations  are  clearly  not  energetically  efficient 


(  Begin  ) 


Fig.  2.  Global  energy  planning  with  varying  mass. 


[37—40],  Hence,  the  method  proposed  in  Ref.  [29],  which  is  based 
upon  the  vehicle  longitudinal  motion  (see  Equation  (1)),  will  be 
used  to  derive  the  most  energy-efficient  speed  profile  for  a  typical 
vehicle  whilst  respecting  the  trip  duration.  However,  this  method 
depends  on  M.  So  the  global  energy  planning  is  formulated  as 
follows: 

Given  a  road-trip  destination,  the  corresponding  road-segment 
speed  limit  profiles,  the  position  of  traffic  light  and  the  road  grade, 
propose  an  energy-efficient  planning  by  considering  that  the 
vehicle  mass  may  vary  during  the  trip. 

To  solve  this  problem,  we  use  the  upper  layer  (global  energy 
planning  layer)  diagram  represented  in  Fig.  2  which  will  be  deeply 
discussed  in  next  sections.  At  the  beginning  of  the  trip,  the  plan¬ 
ning  layer  retrieves  the  destination  provided  directly  by  the  driver, 
the  default  vehicle  mass  which  is  considered  as  the  initial  mass 
and  the  other  vehicle  physical  parameters.  With  the  provided  trip 
destination,  each  road-segment  speed  limits  and  traffic  light  (or 
stop  signs)  position  are  also  retrieved  from  a  GIS.  In  addition,  the 
road  grade  profile  is  obtained  from  this  system  and  with  this  in¬ 
formation,  an  energy-efficient  driving  cycle  is  computed  over  the 
entire  trip  route.  Using  the  computed  driving  profile,  the  electric 
power,  the  energy  consumptions  and  the  anticipated  end-of-trip 
DoD  are  deduced  and  provided  to  the  lower  layer  which  role  is 
to  split  the  driver  power  demand  between  the  battery  pack  and  the 
fuel  cell. 


4.2.  Energy  efficient  driving  cycle 
4.2.1.  Problem  formulation 

To  find  an  energy  efficient  driving  cycle,  we  assume  that  the 
vehicle  with  bounded  acceleration/deceleration  (real  vehicle)  is 
following  a  vehicle  with  no  limitation  (virtual  vehicle).  Since  the 
virtual  vehicle  has  no  limitation,  it  can  follow  precisely  the  speed 
limit  profile  (vi(k),  k  —  1, 2, 3, ....  N}  retrieved  from  the  GIS.  N  is  the 
number  of  samples  over  the  route-trip.  The  states  of  the  real  and 
virtual  vehicles  are  designated  by  X(k)  =  [x(k),v(k)]T  and 
Xv(k)  =  [xv(k),vv(k)]T,  respectively.  The  dynamic  state  equations  of 
the  two  vehicles  are  represented  by  Equations  (9)  and  (12), 
respectively. 


rx(fc+i)i  _  r  x(k)+V(k)AT  i 
[v(fc+l)J  [v(k)  +  u(k)AT-a(k)AT\ 

where 

a(k)  =  }  paACdv2(k)  +gs\n(6)  +  gg  cos(6) 


(9) 


(10) 


u(k)  = 


fin(fc) 

M(k) 


(11) 


AT  is  the  sampling  period.  M(k)  is  the  estimated  vehicle  mass  at 
timestamp  k.  The  algorithm  for  mass  estimation  is  presented  in 
Section  4.4. 


[xv(k  +  l)l  _  rxv(k)+%(k)ATl 

k(k  +  l)J  ~  [  Mk)  \  (12) 

Actually,  (vi(k),  k  =  1,  2, 3, ....  N}  is  the  worst  speed  profile  from 
an  energetic  point  of  view,  as  it  implies  large  accelerations  and 
decelerations  at  the  speed  limit  boundaries.  The  problem  of  finding 
an  energy-efficient  driving  cycle  becomes  one  of  selecting  a  solu¬ 
tion  of  an  optimal  trajectory  tracking  problem  between  the  vehicle 
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under  study  (real  vehicle  with  physical  limitations)  and  the  virtual 
vehicle  by  minimizing  the  following  cost  function: 


Pm  and  the  fuel  cell  power  contribution.  Knowing  the  initial  energy 
Eb(0),  the  optimal  battery  pack  energy  profile  is  given  by: 


Jhl(O)  =  Jhl(N) +  £{/«.(*)} 
k  =  0 

where 

Jhl(N)  =  J[X(JV)  -Xv(N)]rQ(N)[X(JV)  -XV(N)] 


(13)  E*b(k+  1)  =  Eb( 0)  -  JZ([VhVm]^Pm(i)  -  Vb'Pfc)^ 

-IVbVmr^ik+l  )AT 

For  a  given  timestamp  k,  the  theoretical  battery  pack  energy  at 

(14)  the  end  of  the  trip  is  determined  using  the  following  equation: 


and  where 

jHi(k)  =  ~[X(fc)  -Xv(k)]TQ(k)[X(k)  -  Xv(k)]  + 1  u2(k)  (15) 

and  where: 

xv(k+l)  —  xv(k)  +  vv(k)AT  is  the  trajectory  obtained  when  the 
virtual  vehicle  is  moving  along  the  trip  route  at  the  speed  limits  v\. 
Q_  >  0  represents  the  penalty  matrix  associated  with  the  deviation 
between  Xv  and  X. 

The  first  term  of  the  cost  function  (15)  will  let  the  real  vehicle 
trajectory  be  close  to  the  virtual  one  over  the  trip  route,  whereas 
the  second  term  is  introduced  to  take  into  account  the  requirement 
of  low  accelerations/decelerations  on  the  real  vehicle  when 
tracking  the  virtual  vehicle  through  the  speed  limits.  These  low 
accelerations/decelerations  are  one  of  the  key  characteristics  of  an 
energy-efficient  driving  cycle. 

4.2.2.  Optimal  state  tracking  solution 

The  optimal  command  sequence  { u*(k ),  k  =  0,...,  N  -  1}  is 
defined  as  [29]: 

u*  =  arg  minJHL(O)  (16) 

with  the  dynamic  constraint  represented  by  Equation  (9)  and  with 
the  following  additional  nonlinear  constraints: 


vmm  <  <k)  <  vm3X 


where  um in,  Umax,  vmjn  and  vmax  represent  the  allowed  real  vehicle 
minimum  and  maximum  commands,  minimum  and  maximum 
speed  values,  respectively. 

In  the  present  paper,  the  dynamic  programming  method  has 
been  used  to  compute  the  solution  u*(k)  [  10,41 -43|.  The  optimal 
state  sequence  {X*(k)  =  [x*(k),v*(fc)]T,  k  =  0,...,  N]  is  given  by: 


x*(k  + 1)1  _  [  x*(k)+i/*(k)AT  I 

v*(k+  1)  J  [v*(k)  +u*(k)hT  -  a*(/<)ArJ 


(19) 


where  a*(k)  =  (l/2M(k))p3ACd(v*(k))2  +g  sin(fi)  +  pg  cos(6 ) 

The  energy-efficient  driving  cycle  is  therefore  given  by  the 
expression  of  (v*(k),  k  =  1,...,  N}  from  Equation  (19). 


4.3.  Energy-efficient  power  consumption  planning 

Using  the  energy-efficient  driving  cycle  v*  and  the  optimal 
command  u*,  the  corresponding  energy-efficient  power  is  given  by: 

P^(k)  =  M(k)u*(k)v*(k)  (20) 


E£(N)  =  Eb(0)  -  (0  -  VPfc)Ar 


(22) 


It  is  worth  mentioning  that  Eb(N)  is  directly  related  to  the  antici¬ 
pated  battery  pack  depth-of-discharge.  Since  it  is  desirable  to  not  over¬ 
discharge  the  battery  pack,  a  minimum  energy  E™"  should  be  kept.  In 
other  words,  £b(N)  >  Eg1111.  The  fuel  cell  must  be  operated  to  avoid 
violating  this  condition.  One  of  the  lower  layer  goals  is  the  optimal 
operation  of  the  fuel  cell  in  order  to  fulfill  to  the  battery  pack  minimum 
energy  requirement.  This  layer  detail  design  is  present  in  section  5. 


4.4.  Online  vehicle  mass  estimation 

To  be  able  to  predict  the  energy  consumption  accurately,  the 
vehicle  mass  needs  to  be  estimated.  If  we  assume  that  the  vehicle 
dynamics  is  dominated  by  its  longitudinal  motion  and  the  road  grade 
is  known,  the  RLS  method  proposed  in  Ref.  [  23  ]  is  a  good  candidate  for 
the  vehicle  mass  estimation.  This  method  can  be  implemented  in  real¬ 
time  and  is  simple  to  use  in  practice.  Indeed,  the  propulsion  torque  rm 
required  during  the  estimation  process  is  obtained  knowing  the 
electric  power  to  mechanical  power  transformation  map  ijm: 

*m(k)  =  Ptn^"L  (23) 

where  Rm  is  the  wheel  radius  divided  by  the  total  gear  ratio  and 
Pm(k)  is  the  mechanical  power  which  can  be  written  as  a  first  order 
autoregressive  model  with  external  input  (ARX)  and  which  is 
deduced,  when  the  vehicle  speed  v(k)  at  timestamp  k  is  different  to 
zero,  from  Equation  (4)  as  follows:  [29] 

Pm(k)  =  Vm(Pe(kfPm(k-\))  (24) 

As  proposed  in  Ref.  [23]  and  as  we  assume  that  the  trip  route  is 
known  in  advance,  the  route  grade  (0(k),  k  =  1,  2,...,  N]  can  be  also 
known.  Therefore,  the  dynamic  Equation  (1 )  is  rewritten  as  follows: 

Mv(k)  =  -l2PaCdAi?(k) -Mgsin(0(k)) 

-Mgpcos(d(k))  (25) 

where  Jm,  (om  represent  the  powertrain  inertia,  the  powertrain 

rotation  acceleration,  respectively.  Jm  and  Rm  are  provided  by  the 
vehicle  manufacture. 

From  Equation  (25),  the  vehicle  net  acceleration  v  is: 

Kk)  =  ~  ((Tm  ~  \PaCdAAk))  -gsin{d{k)) 

-  gp  cos (6(k))  (26) 


In  practice,  the  driver  may  not  follow  exactly  the  optimal  driving 
cycle.  Hence,  we  must  take  into  account  the  real  mechanical  power 


Rewriting  Equation  (26)  yields  the  following  linear  para¬ 
metric  expression  in  which  the  unknown  parameter  M  can  be 
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estimated  using  the  well-known  recursive  least  squares  (RLS) 
algorithm: 

0(k)  =  mm  (27) 


P}c(k) 


where  the  observation  sequence  O(k),  k  —  1,  2,...,  N  ( N  is  the 
number  of  samples)  is  represented  by: 

0(k)  =  v(k)  +  g  sin (0(k))  +  gp  cos (d(k))  (28) 

and  where  the  regressor  <j>(k)  is  represented  by: 

m  =  ((Tm^<Jm)-^PnQAv2(k))  (29) 

and  where  the  parameter  to  be  identified  is: 

m  =  i  (so) 

The  estimated  mass  M  is  deduced  from  the  recursive  estimation 


pmax 


jjjmax 


-Pb(N) 


5.  Cost  effective  power  splitting 


Fig.  3.  Power  hysteresis  function. 


In  this  section,  we  present  the  detail  design  of  the  lower  layer 
which  is  responsible  for  the  power  demand  splitting  between  the 
fuel  cell  and  the  battery  pack.  At  each  timestamp  during  the  trip, 
the  upper  layer  provides  the  anticipated  final  energy  £f,(N )  in  the 
battery  pack  (see  Equation  (22)).  This  value  is  used  to  compute  the 
DoD.  The  first  requirement  for  the  lower  layer  is  to  keep  the  final 
energy  Eb(N)  greater  than  Eg™.  In  addition,  it  should  take  into  ac¬ 
count  some  desirable  criteria: 

1.  Always  operate  the  fuel  cell  at  its  maximum  efficiency  in  order 
to  minimize  hydrogen  consumption. 

2.  Energy  cost:  since  the  energy  cost  in  battery  may  be  different 
with  hydrogen  cost,  the  splitting  method  must  reduce  the  total 
energy  usage  cost. 

3.  Driving  power  compliance:  the  splitting  method  must  let  the 
vehicle  mechanical  power  to  be  close  to  the  power  demand  in 
order  to  give  to  the  driver,  the  feeling  that  he  is  actively  con¬ 
trolling  the  vehicle  dynamics. 

One  way  to  integrate  these  criteria  when  designing  the  lower 
layer  power  splitting  method  is  to  formulate  and  solve  an  optimal 
control  problem. 

5.1.  Optimal  power  splitting  problem  formulation 

For  each  timestamp  k,  let’s  define  a  fuel  cell  power  reference 
profile  Pfc(k)  based  on  power  hysteresis  function  as  shown  in  Fig.  3: 

•  Pfc(k)  =  Pg?ax  if  E£(N)  <  Eg11" 

.  Prfc(k)  =  OifEg(N)  >Eglax 
.  pr.(k)  =  Pfc(k  -  1)  if  Eg™  <  E*(N)  <  Eglax 

where  Pjgax  is  the  fuel  cell  power  production  at  maximum 
efficiency. 

Minimizing  the  following  cost  function  could  provide  a  good 
trade-off  for  the  splitting  method: 

kL(k)  =i{(Pfc(k)-Pfc(k))2+Cp2(k)  +  (Pm(k)-Pu(k))2}  (31) 

where  C  represents  the  cost  ratio  defined  as  the  ratio  between 
hydrogen  and  grid  electricity  costs. 


The  first  term  of  Ju.  is  added  in  order  to  take  into  account  the 
minimum  battery  energy  requirement.  Indeed,  if  the  upper  layer 
anticipates  that  Eb(N)  will  be  less  than  Eg™,  PL  will  be  equal  to  Pfgax 
and  the  only  way  to  minimize  (Pfc(k)  -  Pfc(k)r  is  set  the  fuel  cell  to 
operate  close  to  its  maximum  power.  The  second  term  takes  into 
account  the  energy  cost  ratio  and  will  be  minimized  if  PfC  =  0. 
Finally,  the  third  term  is  used  for  power  compliance.  Indeed,  if  the 
powertrain  mechanical  power  Pm  is  closed  to  the  driver  power 
demand  Pu,  this  term  is  minimized. 

5.2.  Real-time  power  splitting  solution 

To  allow  the  design  of  a  simple,  practical  and  real-time  solution, 
we  propose  the  following  power  splitting  algorithm  which  pro¬ 
vides  the  local  optimal  control  signals  for  the  fuel  cell  PfC  and  the 
powertrain  P*, : 

•  set  Pg,  to  Pu:  therefore,  the  third  term  is  minimized. 

•  since  the  fuel  cell  power  Pfc  can  take  only  two  different  values  (0 
and  P(gax),  use  each  of  them  to  evaluate  Ju. 

•  select  the  fuel  cell  power  that  gives  the  minimum  value  of  Ju.- 

•  set  Pf  to  that  value. 

6.  Validation  using  experimental  data 

Since  the  energy  planning  layer  is  independent  on  the  way  the 
power  is  splitted  by  the  lower  layer,  we  propose  to  validate  the 
following  aspects  pertaining  to  this  layer: 

•  the  offline  optimal  driving  cycle  calculation  using  dynamic 
programming  and  GIS  data; 


Table  1 

Vehicle  parameters. 

Parameter  Value  Parameter  Value 

M  986  kg  (default  value)  A  2.3  m2 

/)„  1 .22  kg  m  3  fi  0.018 

Cd  0.7  6.  Vk  0  rad 

Umax  4  rn  S  2  Umin  -2  m  S-2 
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Fig.  4.  Driving  cycles  comparison. 


Fig.  5.  Driving  cycles  and  mass:  (a)  experimental  driving  cycle  (blue  curve)  and  optimal  driving  cycle  using  road  segment  speed  limits  (red  curve),  (b)  True  mass  (red  curve)  and 
estimated  mass  using  RLS  algorithm  (blue  curve).  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 
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•  the  online  vehicle  mass  estimation  using  the  consumed  electric 

power,  the  produced  mechanical  power  and  the  RLS. 

In  addition,  the  advantages  of  adapting  energy  planning  method 
with  the  mass  estimation  is  shown  through  a  comparative  study. 

6.1.  Experimental  scenario 

A  small  truck  (maximum  speed  40  Km  h”1)  equipped  with  a 
lead-acid  battery  pack  (rated  to  72  V,  180  Ah)  and  a  fuel  cell  system 
(1.5  lcW)  is  used  for  experiments.  The  sampling  frequency  is  set  to 
10  Hz.  The  vehicle  physical  parameters  are  summarized  in  Table  1. 
One  can  notice  that  the  vehicle  mass  provided  by  the  manufacturer 
is  986  kg  (default  value).  However,  in  our  test  scenario,  the  truck 
mass  at  the  beginning  of  the  experiment  is  measured,  and  it  was 
found  that  the  value  is  1250  kg. 

During  the  experiment,  the  truck  moved  on  a  flat  road  (with  a 
0  grade  value)  and  the  driver  was  instructed  not  to  exceed  the 
speed  limit  of  20  km  h-1  (5.56  m  s-1).  Between  k  =  51  s  and 
k  =  56  s,  the  truck  is  stopped  and  a  load  of  150  kg  added.  Then,  its 
started  moving  again  on  this  flat  road  upto  k  —  117  s. 

6.2.  Optimal  driving  cycle  validation 

Since  the  dynamic  programming  is  a  time-consuming  process, 
the  first  task  is  to  evaluate  the  impact  of  the  varying  mass  on  the 
optimal  driving  cycle.  So  two  driving  cycles  were  computed:  the 
vehicle  default  mass  is  used  for  the  first  one  whereas  a  1400  kg  is 
used  for  the  second  one. 

Fig.  4  shows  the  speed  limit,  the  driving  cycle  if  the  vehicle  mass 
is  986  kg  (default  mass  value)  and  the  driving  cycle  obtained  when 
the  vehicle  mass  is  1400  kg.  The  speed  limit  profile  (dashed  curve) 
has  sharp  value  changes  at  k  =  51  s  and  k  =  56  s.  The  optimal 
driving  cycles  obtained  with  each  of  the  two  masses  are  similar  (the 


blue  curve  (in  the  web  version)  is  similar  to  the  red  one).  The 
relatively  small  root  mean  squared  error  (0.0015)  indicates  that  the 
two  curves  are  very  close  to  each  other.  This  result  suggests  that  at 
low  speed,  the  optimal  driving  cycle  is  little  sensitive  to  the  varying 
mass.  Therefore,  using  the  default  mass  value,  the  optimal  driving 
cycle  is  computed  offline  once.  On  the  other  hand,  the  energy 
consumption  is  dynamically  updated  with  the  vehicle  mass  esti¬ 
mation.  At  high  speed,  the  optimal  driving  cycle  may  be  sensitive  to 
the  vehicle  total  mass.  In  this,  the  optimal  driving  cycle  could  be 
computed  offline  (with  dynamic  programming)  for  the  worst-case 
scenario  using  the  maximum  allowed  mass. 


6.3.  Online  mass  estimation  validation 

Three  different  mass  values  are  involved  in  the  experiment:  the 
default  mass  (986  kg),  the  mass  at  the  beginning  of  the  experiment 
(1250  kg)  and  the  mass  after  k=  56  s  (1400  kg).  The  high  level  layer 
knows  only  the  default  mass  value  before  the  vehicle  started 
moving.  The  vehicle  speed  v  and  the  electric  power  consumption  Pe 
are  monitored  at  each  timestamp  k.  The  power  transformation  map 
77m (k)  —  Pm(fc)/Pe(k)  has  been  estimated  using  the  following 
empirical  expression  [29]: 

Vm(k)  =  0.57Pm(k  -  l)P“1(k)  +  0.2176,  Pe(k)*0  (32) 

Using  the  estimation  presented  in  section  4.4,  the  vehicle  mass 
is  iteratively  evaluated  as  followed: 

•  With  the  measured  value  of  Pe(k)  and  using  Equation  (24),  the 
mechanical  power  Pm(k)  is  estimated. 

•  The  estimated  Pm(k)  and  the  speed  v(k)  is  used  to  deduce  the 
torque  im(k)  (see  Equation  (23)). 

•  Knowing  Jm  and  wm,  the  RLS  algorithm  is  applied  to  find  the 
estimation  M  of  the  vehicle  mass  M. 


4000 

3  2000 


(a) 


Fig.  6.  Power  and  energy  profiles:  (a)  power  profile;  (b)  energy  profile. 
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Fig.  7.  Power  splitting  with  the  default  vehicle  mass. 
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Fig.  8.  Power  splitting  with  adaptive  vehicle 
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Fig.  5,  shows  the  mass  estimation  results.  The  vehicle  speed 
(blue  curve)  and  the  optimal  driving  cycle  (red  curve)  are  repre¬ 
sented  on  the  graph  (a)  of  this  figure.  On  the  graph  (b)  is  shown  the 
estimated  mass  during  the  vehicle  motion  and  the  true  mass. 
Although  the  default  mass  of  986  kg  is  different  with  the  initial 
mass,  one  can  observe  that  M  converges  quickly  towards  1250  kg 
which  is  the  true  mass  for  k  <  56  s.  After  the  mass  has  been  changed 
( k  >  56  s),  the  estimation  algorithm  converges  towards  1400  kg 
which  is  the  true  mass.  Thus,  the  proposed  iterative  estimation 
method  provides  satisfactory  results. 


6.4.  Effect  of  adapting  energy  planning  with  estimated  mass 

As  shown  in  Fig.  2,  we  use  the  estimated  mass  to  dynamically 
update  the  planned  power  usage  as  well  as  the  corresponding  en¬ 
ergy  consumption  and  the  end-of-trip  DoD.  These  two  profiles  are 
presented  in  Fig.  6  where  we  compared  the  power  and  the  energy 
profiles  by  considering  the  vehicle  default  and  estimated  mass.  On 
graph  (a),  we  observe  that  there  is  a  significant  difference  between 
the  power  obtained  using  the  default  mass  value  (blue  curve  (in  the 
web  version))  and  the  one  with  the  estimates  mass  (red  curve).  This 
difference  is  the  main  reason  of  the  under  estimated  energy  con¬ 
sumption  with  the  default  mass  value  (blue  curve  of  graph  (b)).  In 
addition,  we  showed  on  graph  (b)  that  estimating  the  mass  only 
once  at  the  beginning  of  the  trip  is  not  enough  to  predict  the  vehicle 
energy  consumption.  Indeed,  in  our  test  scenario,  the  initial  mass  is 
1250  kg  and  the  dashed  line  on  graph  (b)  indicates  the  energy 
profile  by  considering  this  mass  for  the  whole  trip.  Clearly,  a  sig¬ 
nificant  difference  exists  between  that  energy  profile  and  the 
profile  with  the  estimated  mass.  This  result  indicates  that  a  good 
mass  estimation  is  required  to  properly  plan  the  vehicle  energy 
consumption  for  a  given  trip. 


7.  Power  splitting  method  validation 

This  section  aims  at  providing  evidence  that  using  adaptive 
mass  estimation  improves  the  hybrid  vehicle  energy  management 
system  and  the  associated  power  splitting  method.  We  used  the 
scenarios  described  in  Section  6.1  to  validate  the  power  splitting 
method  and  the  whole  energy  management  system.  The  following 
parameters  are  used:  P™“  =  1500  W,  £b(0)  =  100  Wh, 
Eg™  =  40  Wh,  E™x  =  50  Wh. 

Two  tests  are  performed:  in  the  first  test,  we  compared  the 
battery  pack  (DoD)  when  the  proposed  power  splitting  method  is 
used  with  and  without  adaptive  mass  estimation.  The  DoD  is 
defined  as: 

DoD(k)  =  £b(Ej0)b(k)  (33) 

The  second  test  is  related  to  the  analysis  of  the  impact  of  the 
hydrogen/electricity  cost  ratio  on  the  battery  DoD  with  adaptive 
mass  estimation. 

7.1.  Mass  estimation  impact  on  power  splitting 

The  cost  ratio  C  is  set  to  1  meaning  that  the  hydrogen  and 
grid  electricity  costs  are  similar.  In  addition,  the  fuel  cell  is 
supposed  to  be  stopped  at  the  beginning  of  the  trip.  Fig.  7 
shows  the  result  when  the  default  mass  (986  kg)  is  used  with 
the  power  splitting  method.  On  graph  (a),  the  estimated  battery 
end-of-trip  energy  at  the  (££(N))  using  the  default  known  mass 
(red  curve)  becomes  less  than  Eg110  =  40  Wh  at  k  >  45.7  s.  So 
the  power  reference  profile  P[c(k)  changed  from  0  to  P™ax  and 
the  fuel  cell  was  started  until  the  end  of  the  trip  (see  graph 
(c)).  Although  the  fuel  cell  contributes  to  provide  energy, 
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E*(N )  =  38.71  Wh  which  is  less  than  £™n  =  40  Wh.  The  end- 
of-trip  DoD  is  estimated  to  61.29%. 

Using  the  adaptive  mass  estimation  during  the  trip  (1250  kg  and 
1400  kg),  Fig.  8  shows  that  the  estimated  battery  energy  at  the  end 
of  the  trip  (E£(N))  has  reached  the  minimum  value  (£™n)  in  two 
times  and  so  the  fuel  cell  was  started  twice:  from  k  —  1  to  k  =  23.2  s 
and  from  k  —  50.9  s  to  k=  117.4  s.  The  estimated  end-of-trip  battery 
energy  is  45.24  Wh  which  is  greater  than  Eg™.  Hence  the  minimum 
battery  energy  is  not  violated  with  the  adaptive  splitting  method. 
Furthermore,  the  corresponding  DoD  value  is  54.76%  which  is 
greater  than  the  DoD  value  obtained  with  the  default  mass.  This 
DoD  improvement  has  a  positive  impact  on  the  battery  pack 
degradation  and  efficiency. 


7.2.  Fuel  cost  impact  on  power  splitting 

We  consider  that  the  vehicle  mass  is  continuously  estimated.  To 
analyze  the  impact  of  the  energy  cost  ratio  on  the  DoD,  we  run  several 
tests  in  which  C  value  varies  from  0.1  (hydrogen  is  much  cheaper  than 
electricity)  to  4  (hydrogen  is  much  more  expensive  than  electricity). 
Fig.  9  shows  the  DoD  as  a  function  of  C.  One  can  observe  that  as  the  cost 
ratio  is  less  or  equal  to  1  (hydrogen  is  less  or  as  expensive  as  elec¬ 
tricity),  the  splitting  method  allows  the  fuel  cell  to  contribute  to 
vehicle  traction  power.  Indeed,  from  Equation  (32),  by  setting  Pm  =  Pu, 
Julk)  becomes:  JLL(k)  =  l/2{(Pfc(k)  -  Pfc(k))2  +  CPfc(k)}.  The  rela¬ 
tive  weight  ofp2  is  less  than  the  weight  of  (Pfc(k)  -  P'fc(k))2.  One  way 
to  minimize Jll  is  to  let  PfC  to  follow  Pfc.  On  the  other  hand,  if  Cis  greater 
than  1  meaning  that  hydrogen  is  more  expensive  than  electricity,  the 
relative  weight  of  (Pfc(k)  -  P[c(k))2  is  less  than  the  weight  of  Pfc.  So  to 
minimize Jn,  the  splitting  method  sets  Pfc  to  0.  Therefore,  no  hydrogen 
is  used. 


8.  Conclusion 

We  proposed  an  adaptive  energy  management  and  power 
splitting  method  for  fuel  cell  hybrid  electric  vehicle.  This  method  is 
based  on  a  two  layer  architecture  in  which  the  upper  one  is 
responsible  for  estimating  the  battery  pack  depth-of-discharge 
using  the  online  mass  evaluation.  Given  the  estimated  depth-of- 
discharge,  the  lower  layer  uses  an  optimal  power  splitting  algo¬ 
rithm  to  share  the  power  demand  between  the  fuel  cell  and  the 
battery  pack.  The  goal  of  this  splitting  algorithm  is  to  improve  the 
depth-of-discharge  while  taking  into  account  the  cost  ratio 
hydrogen/electricity.  In  addition,  the  fuel  cell  is  powered  on  to  al¬ 
ways  work  at  its  maximum  efficiency.  A  comparative  study 
demonstrated  that  by  taking  into  account  the  online  mass  esti¬ 
mation,  the  overall  energetic  efficiency  is  improved. 
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